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ABSTRACT 

Finite-amplitude  standing  waves  in  air  at  ambient  con- 
ditions contained  in  a  rigid-walled  cylindrical  tube  with 
a  large  length-to-diameter  ratio  were  experimentally  inves- 
tigated.  The  pressure  waveform  at  the  end  of  the  tube  was 
digitized  and  Fourier  analyzed  on  an  IBM  360  digital 
computer.   Amplitudes  and  phases  were  obtained  for  all 
harmonics  with  amplitudes  greater  than  1%  of  the  fundamental 
for  strength  parameters  from  0.25  to  1.00  and  for  frequency 
parameters  from  -0.8  to  2.0.   The  strength  and  frequency 
parameters  are  defined  as  r4bQ  and  2Af/Afj  respectively, 

'2 

where  M  us  tne  Mach  numtjer  Oo.  funv.iamental  ,  i->   the  ncnj-iLn~ 
iarity  parameter,  Q  the  quality  factor  of  the  resonator, 
A  f  the  frequency  away  from  fundamental  resonance ,  and 
A  f  1  the  band  width  at  the  half-pov;er  points.   V^hen  these 
results  are   compared  to  the  theoretical  model  of  Coppens 
and  Sanders,  it  is   seen  that  v;hile  the  theory  accurately 
predicts   the  magnitude  and  shape  of  the  harmonic  content, 
it  consistently  underestimates  the  frequency  at  which  each 
harmonic  peaks.   In  addition,  the  theory  fails  to  predict, 
except  qualitatively,  the  phase  angles. 


TABLE  OF  CONTENTS 

I.  INTRODUCTION  5 

II.  EXPERIMENTAL  CONSIDERATIONS ■ 9 

A.  APPARATUS 9 

B.  MICROPHONE  CALIBRATION  10 

•  C.   MICROPHONE  OUTPUT  AND  STRENGTH  PARAMETER  12 

D.  ATTENUATION  CONSTANT ■ 12 

E.  EREQUENCY  PARAMETER  13 

III.  COMPUTER  ANALYSIS 14 

A.  SAMPLE  RATE 14 

B.  THE  DIGITAL  COMPUTER  PROGRAM 15 

C.  THE  ANALOG  TO  DIGITAL  CONVERSION 16 

D.  SYSTEM  TEST 18 

IV.  PROCEDURES 19 

A.  PRERUN  PROCEDURES ■ 19 

B.  DATA  TAKING  PROCEDURES 20 

V.  RESULTS 21 

APPENDIX  A:   A-D  PROGRAI^  PARAJ^.ETFRS  AND  OPTIONS 2  3 

APPENDIX  B:   EIGURES  AND  PLOTS 26 

COMPUTER  PROGRAMS • 39 

BIBLIOGRAPHY ■ ^^ 

INITIAL  DISTRIBUTION  LIST  ^^ 


FORM  DD  1473 


52 


LIST  OF  FIGURES 

1.  Apparatus  Block  Diagram  26 

2.  PI  6200  Dynamic  Range 27 

3.  Analog  Patch  Board ■ 28 

4.  Logic  Patch  Board  29 

5.  Test  Signal  Analysis 30 

6.  Point  Scatter ^1 

7.  Harmonic  Amplitudes  SP  =  0.25 —  32 

8.  Harmonic  Amplitudes  SP  =  0.50 33 

9.  Phase  Plot  SP  =  0.50 34 

10.  Harmonic  Amplitudes  SP  =  0.75 35 

11.  Phase  Plot  SP  =  0.75  36 

12.  Harmonic  Amplitudes  SP  =  1.00 37 

13.  Phase  Plot  SP  =  1.00  


38 


I.   INTRODUCTION 

Coppens  and  Sanders  [1/2]  have  developed  a  one- 
dimensional,  nonlinear  wave  equation  with  dissipative  terms 
corresponding  to  those  encountered  in  a  rigid-walled  closed 
tube  with  a  large  length-to-diameter  ratio.   An  attenuation 
constant  of  the  form 


n     1 


is  assumed,  where  a   is  the  absorption  coefficient  associa- 

n  ^ 

ted  with  the  frequency  nw  which  characterizes  the  v/all 
losses,  R  is  the  radius  of  the  tube  and  L  is  the  length  of 
the  tube.   With  this  additional  loss  mechanism  the  v/ave 
equation  takes  the  form 
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Y  =  C  /C   the  ratio  of  specific  heats, 
'      p   v  ^ 

V  =  kinematic  viscosity, 
Pr  =  Prandtl  number, 

w   =  angular  frequency  of  the  fundamental, 
n   =  harmonic  number, 


C   =  the  thermodynamic  speed  of  sound, 

a-j^  =  absorption  coefficient  for  the  lowest  frequency 

and 

U  =  particle  velocity  =  ttt-  +  S  U  . 
^  -^    9t      n 

n 

To  arrive  at  the  above  equation  it  was  necessary  to  make  the 
simplifying  assumptions  that  the  wall  losses  were  much 
greater  than  the  bulk  losses  and  that  the  Mach  number 

M  =  P,/p  C  ^ 
1   o  o 

is  much  less  than  unity.   P,  is  the  peak  value  of  the  funda- 
mental of  the  pressure  at  the  rigid  end  of  the  tube. 
A  Fourier  series  solution  of  the  form 

P(x.t)  CO  ^         ^^ 

~-~  ^   M      S      R(n)    Co£|n7T(l-    ^)    Sin  (nfot  +  ffm ) 

p  C        n=l 
o  o 

is  assumed  with  the  tube  excited  with  a  pressure  distribu- 
tion given  by 

P,  (x,t)  =  Mp  C  ^  CosTT  (l-x/L)Sincot. 
11       ^o  o 

When  this  solution  was  substituted  directly  into  the  non- 
linear wave  equation,  the  following  set  of  algebraic 
equations  resulted: 

and 
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infinitesimal  amplitude  absorption  coefficient 
in  free  space, 

angular  frequency  at  resonance, 

strength  parameter 


Mbk/2  T  =  PTbk/(2p  C   a,) 
1     1       o  o   1 


SP 

k   =  wave  number 


and 


Aw  =  to  -  to  . 
r 


These  equations  were  then  programmed  for  computer  solution 

For  given  values  of  Ato  and  SP ,  the  computer  calculated 

values  of  P  /P,  and  Q    .      At  present,  solutions  have  been 
n   1       n      ^ 


found  for  strength  parameters  up  to  0.75. 


Previous  investigations  in  this  area  have  been  carried 
out  by  Sanders  [1,2]  and  Beech  [3].   Both  investigators 
examined  the  harmonic  content  of  the  pressure  waveform  at 
the  rigid  end  of  the  tube  near  resonance  conditions. 
Sanders  employed  a  wave  analyzer  to  obtain  the  frequency 
response  for  the  first  four  overtones.   Good  agreement  with 
theory  was  found  but  no  phase  information  was  obtained. 
Beech  used  photographic  and  graphical  techniques  to  provide 
numerical  analysis  of  pressure  waveforms  near  resonance. 
Phase  information,  however,  was  obtained  only  for  cases  at 
resonance.   Again  good  agreement  with  theory  was  obtained 
for  cases  approaching  the  onset  of  shock. 

The  present  investigation  examines  both  the  harmonic 
content  and  phase  relationships  for  strength  parameters  as 
high  as  0.75.   Again  numerical  methods  were  employed. 
Whereas  the  graphical  techniques  of  Beech  were  used  to 
generate  64  data  equally  spaced  throughout  the  pressure 
wave  form,  the  technique  in  this  investigation  employed 
analog-to-digital  conversion  by  computer  to  provide  better 
than  1000  data  to  represent  the  pressure  waveform.   With 
the  pressure  waveform  digitized  to  such  anextent,  it  was 
possible  to  obtain  amplitude  and  phase  information  for  all 
harmonics  with  amplitudes  greater  than  one  percent  of  the 
fundamental . 


II.   EXPERIMENTAL  CONSIDERATIONS 

A.   APPARATUS 

A  block  diagram  of  the  experimePxtal  system  used  in  this 
investigation  is  shown  in  Figure  1.   Finite-amplitude 
standing  waves  were  generated  in  air  at  ambient  conditions 
contained  within  a  tube  six-feet  long  with  an  inner  diameter 
of  2.250  in.  and  a  v;all  thickness  of  1.125  in.   The  air  was 
excited  at  one  end  by  a  piston  which  was  driven  by  an  M-B 
Electronics  Model  EA  1500  exciter  capable  of  a  maximum 
no-load  acceleration  of  124g.   The  EA  1500  was  driven  by  two 
M-B  electronics  Model  2120MB  power  amplifiers  operating  in 
parallel.   The  maximum  attainable  acceleration  v;ith  the 
weight  of  the  piston  as  load  was  about  50g.   The  frequency 
and  power  level  were  controlled  by  a  General  Radio  Type 
1161-AGC  coherent  decade  frequency  synthesizer.   With  this 
synthesizer  it  was  possible  to  select  frequency  to  better 
than  0.001  cycles/sec.  (Monitoring  the  frequency  v/ith  a 
Hewlett-Packard  Model  521C  frequency  counter  showed  the 
frequency  stability  to  be  better  than  one  part  in  10  million.) 
The  acceleration  of  the  piston  was  determined  by  measuring 
the  output  an  Endeveco  Model  2215  accelerometer  implanted 
within  the  piston.   When  the  piston  was  v/ithdravm  from  the 
tube  and  driven  in  free  air  at  the  approximate  experimental 
frequency  and  acceleration,  the  harmonic  distortion  was 
measured  to  be  less  than  0.3  percent.   The  other  end  of  the 


tube  was  capped  with  a  thick  plate  which  was  securely  bolted 
to  the  tube.   The  sound  pressure  level  at  the  rigid  end  of 
the  tube  was  detected  with  a  1/4-in.  diameter  Bruel  and 
Kjaer  Type  4136  condenser  microph^one  inounted  in  the  cap  so 
that  its  diaphragm  was  flush  with  the  end  of  the  tube.   The 
microphone  response  was  flat  to  within  +0.5  dB  from  80  Hz 
to  20  kHz.   The  outputs  of  the  microphone  and  accelerometer 
were  monitored  with  Hev/lett-Packard  400D  vacuum-tube  volt- 
meters.  The  microphone  output  also  went  to  a  Hewlett-Packard 
Model  302A  wave  analyzer  so  that  the  first  tV70  overtones 
could  be  monitored.   The  microphone  output  was  recorded  by 
a  PI  6200  General  Purpose  Portable  Instrumentation  Tape 
Recorder.   Recordings  were  made  in  the  FM  mode  at  a  tape 
speed  of  37.5  in/s  so  that  the  frequency ' response  of  the 
tape  recorder  was  +  1  dB  from  DC  to  10  kHz.   Although  the 
signal  was  recorded  at  37.5  in/s  it  was  necessary  when 
digitizing  the  signal  to  play  back  at  a  speed  of  3.75  in/s. 
To  insure  that  this  change  in  tape  speeds  did  not  alter 
the  recorded  frequencies,  a  test  signal  was  recorded  at  the 
higher  speed  and  then  played  back  at  the  slow  speed.   The 
two  signals  were  observed  on  a  dual  trace  oscilliscope 
where  it  was  observed  that  the  frequency  of  the  test  signal 
was  10.0  times  that  of  the  signal  played  back  at  3.75  in. 
per  second  with  no  discernable  drift  between  the  tv70  signals. 

B.   MICROPHONE  CALIBRATION 

Investigation  of  the  standing  waves  at  various  strength 
parameters  requires  an  absolute  calibration  of  the  microphone 
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system.  The  microphone  system  is  defined  here  as  the  micro- 
phone, preamplifier  and  the  equipment  loaded  onto  the  output 
of  the  preamplifier.   The  system  sensitivity  is  defined  as 

S   =  V  /P 
m    m'^ 

where  V   is  the  peak  voltage  output  of  the  preamplifier  and 
P  Is  the  peak  pressure  at  the  microphone.   Tv/o  methods  v;ere 
employed  to  determine  S^.   The  first  and  simplest  used  a 
Bruel  and  Kjaer  Model  4220  pistonphone  which  produces  a 
known  sound  pressure  level  at  250  Hz.   The  sensitivity  v;as 
found  to  be  S   =  (-77.1  +  0.2)dB  re  1  volt/ybar,  where  the 
uncertainty  is  that  inherent  in  the  pistonphone  specifications 
The  second  method  makes  use  of  theoretical  predictions  by 
Coppens  for  small  strenc^th  pa,ra;Be-ters  as  verified  by  Beech 
[3] .   Coppens  shows  that  at  resonance  the  ratio  of  the  peak 
pressure  of  the  second  harmonic  to  the  square  of  the  peak 
pressure  of  the  fundamental  is  given  by 

P2/P^  =  (v/2~  w6)/(8p^C^a)  . 

\Vhen  the  definition  of  sensitivity  is  substituted  into  the 
above  equation  the  result  is 

S^  -  (6a3V^)/(4/2~  p^C^aV2), 

S_  was  then  found  by  exciting  the  tube  at  resonance  and 

measuring  the  fundamental  and  second-harmonic  voltages.   By 

this  method  S   =  (-77.2  +  0.3)dB  re  1  volt/ybar  where  the 
m  — 

uncertainty  is  due  largely  to  the  error  in  determining  the 
attenuation  constant. 
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C.   MICROPHONE  OUTPUT  AND  STRENGTH  P/^RAMFTER 

In  this  investigation,  where  the  pressure  waveform  was 
to  be  analyzed  at  various  strength  parameters ,  it  was 
necessary  to  find  a  relationship  between  the  strength 
parameter  and  the  observable  quantities.   Using  the 
equations 


I  V   ^ 

S   =  SPL  +  21  log  -=^=-F 
m  ^1.775 


+  72.0, 


P       P      P 

M   =  — i^  =  p^  -  — ^-2  =  SPL  -  197.0, 

p  C       o  p  C 
o  o  o  o 


and 


M  =   SP(2a/K) 
b 

one  can  show  that  the  rms  valae  of  the  f undainental  of  the 
microphone  output  (v)  is  related  to  strength  parameter  and 
the  attenuation  constant  by 

V  =  136.0  a  •  SP. 

The  constant  136.0  is  relatively  insensitive  to  small 
temperature  variations  about  room  temperature. 

D.   ATTENUATION  CONSTANT 

The  attenuation  constant  for  infinitesimal  amplitudes 
has  been  shown  by  Beech  to  be 

\ 


IT 

«1  =  c 

l  o  I 


Af, 


where  Af  is  the  frequency  difference  betv;een  the  half  power 
points  on  the  resonance  curve  (taken  for  constant  accelera- 
tion amplitude  of  the  piston) . 
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E.   FREQUFNCY  PARAMETER 

The  frequency  parameter  is  defined  by 


FP  =--  2Aa)/(jo  6^  , 
r  1 


Ao)  =  w  -  0) 


where 


and 

6  ,  =  2  a ,  C  /to  . 
1      lor 

This  frequency  parameter  is  related  to  the  experimental 
frequency  increment  (Af   )  by 

Af   •=  Af  •  FP/2, 
ex 

where  Af  is  the  frequency  difference  between  the  half  power 
points  of  the  infinitesimal  resonance  curve  as  used  in 
calculating  a,  in  this  investigation  Af  was  found  to  be 
1.15  +  O.OIH   so  that 

Af    =  0.58  FP. 
ex 
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III.   COMPUTER  ANALYSIS 

A.   SAMPLE  RATE 

Many  computer  programs  exist  which,  when  supplied  with 
data  points  defining  a  periodic  waveform,  will  give  the 
Fourier  coefficients.   The  first  consideration,  generally, 
is  the  num.ber  of  data  required  per  cycle.   Theoretically 
this  number  should  be  at  least  twice  the  number  of  harmonics 
desired.   However,  the  question  of  accuracy  of  the  computer 
pjredictions  of  harmonic  amplitudes  arises  im.mediately .   As 
might  be  expected,  the  accuracy  of  the  harmonics  determined 
tails  off  faster  as  fewer  samples  per  cycle  are  used.   This 
then  leads  to  the  problem,  of  first  dctcrm.ining  how  m.any 
significant  harmonics  are  desired.   For  finite-amplitude 
standing  v/aves  of  large  strength  parameter  the  waveform 
qualitatively  approaches  a  function  whose  harmonics  content 
falls  off  inversely  with  harmonic  number.   If  such  a  v^ave- 
form  is  passed  through  amplifiers  and  storage  devices  before 
being  analyzed,  harmonics  below  a  given  amplitude  will  be 
lost  in  the  inherent  noise.   The  device  in  this  investiga- 
tion which  was  most  seriously  restrictive  in  this  respect 
was  the  PI  6200  FM  tape  recorder.   It  was  found  that  the 
dynamic  range  of  this  instrument  was  about  42  db  which 
restricted  the  investigation  to  harmonics  with  amplitudes 
greater  than  about  one  percent  of  the  fundamental.   For  a 
periodic  ramp  function  this  restriction  occurs  at  about  the 
30th  harmonic. 
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Determination  of  the  number  of  samples  per  cycle  to  be 
used  was  based  on  the  arbitrary  criterion  that  the  error 
in  the  30th  harmonic  be  less  than  0.1  percent  for  an  ideal 
ramp  functirm.  Applying  several  computer  algorithms  for 
determining  Fourier  coefficients  to  an  ideal  ramp  function 
one  finds  that  1000  samples  per  cycle  are  required  for  the 
above  accuracy. 

B.   THE  DIGITAL-COMPUTER  PROGRAM 

The  primary  function  of  this  program  v;as  to  generate 
normalized  Fourier  coefficients  for  a  periodic  function 
given  digitized  data  on  a  seven-track  tape.   Secondary  to 
this,  the  program  changed  the  form  of  the  coefficients  to 
A(n)  Sin  (Ntot  +  ({)  (n)  and  computed  A(n)/A(l)  for  each  harmonic. 
Optional  functions  provided  by  the  program  were  the  averag- 
ing of  several  periods  of  data  to  reduce  the  affects  of 
random  noise.   An  optional  plot  of  the  input  data  was  avail- 
able along  with  points  reconstructed  from  the  calculated 
Fourier  coefficients.   This  plot  provided  a  visual  check  on 
the  input  data  and  the  computed  coefficients. 

The  digitized  data  provided  by  the  SDS  9300  were  written 
on  seven-track  tape.   However,  the  binary  representation  of 
numbers  on  the  seven-track  tape  differs  from  those  on  the 
standard  nine-track  tape  used  by  the  IBM  360  system.   To 
accomodate  the  seven-track  tape,  a  conversion  program 

(subroutine  FORM)  was  used  to  convert  the  data.   The  con- 
verted data  were  not  floating  point  but  integers  and  required 

further  conversion  v;ith  the  use  of  a  scaling  factor. 
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Prime  consideration  was  given  to  the  types  of  algorithms 
commonly  used  in  generating  Fourier  coefficients.   Two  of 
the  more  widely  used  methods  are  the  Cooley-Tukey  [5]  Fast 
Fourier  transform  and  the  recursive  technique  of  Goertzel 
[4],   Both  methods  were  tested  for  accuracy  in  predicting 
the  coefficients  of  known  waveforms.   The  two  methods  were 
found  to  be  equal.   One  advantage  of  the  Cooey-Tukey  method 
is  its  minimal  use  of  computer  time.   Its  computation  time 
is  about  an  order  of  magnitude  less  than  that  required  by 
the  Goertzel  method.   However,  when  this  time  saving  is 
compared  with  the  overall  required  computer  time  for  analysis 
the  advantage  diminishes.   The  final  consideration  and,  in 
the  end  the  most  important,  was  the  flexibility  in  using 
data.   The  Cooley-Tukey  m.cthod  required  an  integer  pov.'er-of- 
two  number  of  points.   To  be  able  to  satisfy  this  condition 
along  with  the  requirement  for  1000  or  more  samples  per  cycle 
complicated  the  digitizing  process.   The  Goertzel  method 
required  only  an  odd  number  of  samples  and  for  this  reason 
was  used  in  the  computer  program. 

C.   THE  ANALOG  TO  DIGITAL  CONVERSION  (A-D)  • 

The  conversion  of  the  analog  data  was  accomplished  with 
the  use  of  a  hybrid  computer  consisting  of  the  SDS  9300 
high-speed  digital  computer  and  the  Ci  5000  general-purpose 
analog  computer.   Use  of  the  system  for  A-D  conversion 
requires  a  FORTRAN  program  for  the  SDS  9300  and  properly 
constructed  analog  and  logical  patch  boards  for  the  Ci  5000. 
In  addition,  several  peripheral  equipments  are  required. 
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These  are  the  card  reader,  magnetic  tape  transport,  line 
printer  and,  if  desired,  a  pen  recorder. 

Aside  from  the  question  of  sampling,  which  was  deter- 
mined earlier,  the  important  question  in  the  A-D  conversion 
is  accuracy.   The  clock  in  the  computer  is  accurate  to  one 
part  in  100  million,  well  within  any  accuracy  requirements. 
The  A-D  converter  is  designed  to  accept  signals  between 
+  100  volts  and  divide  this  interval  into  16,384  parts 
providing  a  resolution  of  12.2  mV.   For  a  properly  amplified 
analog  signal,  the  dynamic  range  of  the  A-D  converter  is 
much   greater  than  the  PI  6200  tape  recorder. 

The  analog  patch  board  (Figure  3) ,  as  used  for  A-d  con- 
version, served  to  amplify  the  analog  signal  prior  to  con- 
version and  distributed  the  signal  to  the  converter  and 
oscilliscope .   Additionally,  a  patch  was  provided  linking 
the  Digital-to-Analog  converter  with  a  pen  recorder.   This 
additional  function  is  explained  in  detail  in  the  pro- 
cedures section. 

The  logic  patch  board  (Figure  4),  had  three  functions. 
The  first  provided  a  sampling  frequency  (10,000  Hz)  to  the 
A-D  converter.   This  frequency  represents  the  approximate 
maximum  sampling  rate  allowed.   Since  the  approximate 
;i:esoaa,nt  freqriency  of  the  tube  was  94  Hz  it  was  necessary  to 
play  back  signals  from  the  PI  6200  at  one  tenth  the  recorded 
speeds  to  acquire  the  desired  number  of  samples  per  cycle. 
The  second  function  provided  a  triggering  voltage  through 
an  external  switch  (DS-1)  which  was  used  to  start  the  data 
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conversion  process.   Associated  with  this  function  is  an 
externally  controlled  delay  (DFOO)  which  insured  that  the 
triggering  voltage  remained  on  long  enough  for  one  set  of 
data  (one  record)  to  be  written  on  the  seven-track  tape. 
The  last  function  provided  another  triggering  voltage 
through  an  external  switch  (DS-2) .   When  the  last  record  was 
to  be  written  on  tape  switch  DS-2  was  turned  on  while 
switch  DS-2  v/as  triggered. 

D.   SYSTEM  TEST 

Prior  to  analyzing  experimental  data,  it  was  necessary 
to  test  the  entire  system  process  of  A-D  conversion  and 
computer  analysis.   The  results  of  this  test  for  a  ramp 
function  is  shown  in  Figure  5. 


IV.   PROCFDURFS 

A.   PRERUN  PROCEDURES 

Beech  pointed  out  that  the  resonant  frequency  of  the 
tube  changed  by  0.16  Hz  for  each  degree  centigrade  of 
temperature  change.   To  minimize  this  affect,  all  data 
gathering  runs  were  made  during  the  evening  hours  v/hen 
the  room  temperature  remained  relatively  constant  for 
several  hours.   Typical  resonant  frequency  variations  during 
these  hours  were  of  the  order  of  0.01  Hz  per  hour.   Several 
hours  prior  to  taking  data  the  equipment  was  turned  on 
and  allowed  to  stabilize.   During  this  warmup  period,  the 
PI  6  200  tape  recorder  was  calibrated  to  insure  that  the 
FM  system  was  properly  aligned  and  to  insure  that  the 
signals  being  recorded  were  not  overdriving  the  tape  recor- 
der preamplifier  and  thereby  introducing  unwanted  signal 
distortion.   Next,  the  alignment  of  the  piston  in  the  free 
end  of  the  tube  was  checked  by  observing  the  output  of  the 
accelerometer  on  an  oscilloscope  and  measuring,  with  a 
wave  analyzer,  the  second-  and  third-harmonic  content.   If 
necessary,  the  harmonic  content  could  be  reduced  by  carefully 
adjusting  the  alignment  of  the  piston  so  that  the  distortion 
from  any  harmonic  was  less  than  2  percent.   To  minimize 
leakage  of  sound  from  the  tube,  the  piston  was  fitted  with 
an  "O"  ring  and  lubricated.   The  harmonic  distortion  intro- 
duced by  the  piston  was  sensitive  to  the  type  of  lubricant 
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used.   It  was  found  that  a  light  turbine  oil  minimized  the 
distortion  while  providing  an  adequate  seal. 

Just  prior  to  taking  each  set  of  data,  it  was  necessary 
to  determine  the  resonant  frequency  and  the  attenuation 
constant.   For  this  measurement  the  piston  was  excited  at 
a  constant  value  of  acceleration  low  enough  so  that  the 
strength  parameter  at  resonance  was  less  than  0.03.   From 
the  frequency  difference  between  the  half  power  points,  the 
attenuation  constant  was  determined.   The  resonance  frequency 
was  taken  to  lie  half  way  between  the  half  power  points. 

B.   DATA  TAKI-NG  PROCEDURE 

Corresponding  to  the  desired  strength  parameter  was 
a  given  value  of  the  pressure  fundamental  as  shown  earlier. 
With  the  piston  excitation  set  to  the  proper  level  to  keep 
the  pressure  fundamental  constant,  data  v;ere  recorded  at 
0.1  Hz  intervals  for  twenty  frequencies  from  0.6  Hz  below 
resonance  to  1.3  Hz  above.   At  each  frequency,  the  signals 
from  the  microphone  and  accelerometer  were  recorded  on  the 
PI  6200  at  high  speed  (37.5  in/s)  for  about  10  seconds 
(30  ft.),  and  the  second-harmonic  amplitude  of  the  pressure 
was  noted  to  be  later  compared  with  the  computer  analysis. 
Since  the  taking  of  an  entire  set  of  data  at  one  strength 
parameter  required  from  30  minutes  to  an  hour,  at  the  end 
of  the  run  it  was  necessary  to  check  the  resonant  frequency 
and  attenuation  constant  to  insure  that  the  variation  in 
resonant  frequency  v/as  substantially  less  than  one  frequency 
interval  (0.1  Hz)  and  that  the  change  in  attenuation  constant 
was  less  tlian  one  percent. 
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The  theoretical  predictions  (8)  are  compared  with  the 
experimental  results  in  Figures  6  through  13  for  each  value 
of  strength  parameter,  plots  are  shown  of  the  harmonic 
amplitude  vs.  frequency  parameter  and  of  the  phase  difference 
between  consecutive  harmonics  vs.  inverse  harmonic  number. 
Figure  6  shows  the  actual  data  obtained  for  three  different 
experimental  runs.   The  degree  of  scatter  is  a  demonstra- 
tion of  the  excellent  reproducibility  of  the  results.   Both 
experimental  and  theoretical  results  are  limited  to  harmonics 
with  amplitude  greater  than  1%  of  the  fundamental. 

In  comparing  the  experimental  values  of  the  harmonic 

amplitudes  with  the  theoretical  predictions  for  a  strength 

parameter  of  0.25  it  is  seen  that  the  agreement  is  fairly 

good  but  that  there  does  seem,  to  be  a  slight  systemat  c 

difference.   The  corresponding  curves  for  higher  strength 

parameters  display  this  difference  to  a  more  marked  degree. 

It  is  seen  that  the  experimental  results  peak  at  about  the 

predicted  values  of  P ,  x/P/ts  and  that  the  shape  of  the 
^'  (n)'  (1) 

curve  for  each  harmonic  is  about  as  predicted.   However, 
the  experimental  results  peak  at  a  value  of  frequency 
parameter  significantly  greater  than  predicted  and  this 
difference  increases  with  harmonic  number.   It  is  interest- 
ing to  note  that  at  the  frequency  of  maximum  distortion  of 
the  pressure  v7aveform  (at  approximately  FP  =  0.5)  theory 
accurately  predicts  the  harmonic  amplitudes. 

The  comparison  between  the  experimental  and  theoretical 
phase  angles  displays  an  even  greater  discrepancy.   VJhile 


values  of  4)  -    ^      ,  do  seem  to  be  linearly  related  to  1/n , 

n     n-1  ^  '     ' 

the  values  for  the  second  harmonic  are  greatly  different 
and  the  slopes  of  the  experimental  curves  are  significantly 
greater  than  predicted. 

In  conclusion,  it  appears  that  while  the  present  theo- 
retical model  is  able  to  predict  the  magnitude  and  shape  of 
the  harmonic  content  of  the  distorted  waveform,  it  consis- 
tently underestimates  the  frequency  at   which  each  harmonic 
peaks.   This  underestimatior;  .lncif?<j£es  with  harmonic  numl)er 
and  strength  parameter.   In  addition,  it  fails  to  predict, 
except  qualitatively,  the  phase  relations  in  the  distorted 
wave . 
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APPENDIX  A 
THE  A-D  PROGRAM  PARAMETERS  AMD  OPTIONS 

The  purpose  of  this  program  was  to  convert  up  to  four 
channels  of  analog  data  into  a  series  of  digital  points  on 
seven-track  magnetic  tape.   The  program  had  several  input 
parameters  which  had  to  be  provided  to  the  computer  via 
the  teleprinter  after  the  FORTRAN  program  had  been  read  into 
the  system. 

NSAMP  (input  1.)  was  an   integer  which  corresponded  to  the 
number  of  data  or  samples  to  be  taken  from  each  channel. 
For  this  investigation  a  sampling  rate  of  10  kHz  was  used 
and  NSAMP  was  chosen  to  span  one  cycle  of  the  input  waveform, 

NCHAN  (input  2.)  was  an  integer  which  specified  the 
number  of  channels  of  analog  data  being  converted. 

NREC  (input  2.)  specified  hovs'  many  records  were  to  be 
written  on  tape  each  time  sampling  was  initiated.   A  record 
was  defined  as  the  data  set  corresponding  to  the  product  of 
NCHAN  and  NSAMP.   Because  of  the  high  sampling  rate  used  it 
was  necessary  to  set  this  value  to  1. 

ITAPE  (input  4.)  designated  which  magnetic-tape  trans- 
port held  the  seven-track  tape.  The  integer  used  here  had 
to  correspond  that  selected  on  the  "unit  select"  switch  on 
the  particular  tape  unit  used. 

NDEL  (input  5.)  determined  how  rapidly  the  data  on  the 
seven-track  tape  was  read  back  when  using  program  option  7. 
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An  example  of  the  input  parameters  as  typed  on  the 
teleprinter  might  be 

NREC  =  1,  NSAMP  =  1000,  NCHAN  =  3,  ITAPE  =  1, 
NDEL  =  20*   (carriage  return) . 

The  control  symbol  *  indicated  to  the  computer  that  the 
input  parameters  had  been  set.   Care  was  taken  when  choos- 
ing NSAMP  and  NCHAN  to  insure  that  their  product  did  not 
exceed  the  first  dimension  of  IBUF  found  in  the  first  FORTRAN 
program  statement. 

When  the  parameter  list  had  been"  given  the  teleprinter 
responded  with 

"OPTION  =  (II) " . 

The  computer  v/as  then  given  a  one-digit  number  ccrrGspcnding 
to  the  list  of  available  options. 

RESTART  (option  1.)  was  used  whenever  any  one  or  all  of 
the  input  parameters  was  to  be  changed. 

GO  (option  2.)  was  used  when  analog  to  digital  conver- 
sion was  to  be  done.   Upon  receiving  this  option  the  compu- 
ter awaited  the  initiating  trigger  from  switch  DS-1  to 
start  the  actual  conversion  process.   Each  trigger  from  DS-1 
caused  one  record  of  data  to  be  written  on  tape.   When  a 
sufficient  number  of  records  v/ere  written,  DS-1  was 
triggered  with  DS-2  in  the  "up"  position.   This  caused  the 
last  record  to  be  written  followed  by  a  teleprinter  response 

"OPTION  =  (II)  "  . 
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END  OF  FILE  (option  3.)  was  used  to  put  a  mark  on  the 
tape  to  separate  or  isolate  groups  of  records.  Groups  of 
records  separated  by  END  OF  FILE  marks  are  called  files. 

SKIPFILES  (option  5.)  were  used  when  files  already  on 
tape  were  to  be  skipped  over.   When  used,  a  four-digit 
number  was  given  the  computer  to  specify  the  number  of 
files  to  be  skipped. 

LIST  ONE  RECORD  (option  6.)  was  used  in  reading  data 
from  the  seven-track  tape.   Use  of  this  option  caused  data 
written  on  tape  to  be  written  on  the  page  printer  after 
having  specified  the  number  of  data  to  print. 

PLOT  ONE  RECORD  (option  7.),  when  called,  caused  one 
record  to  be  read  from  the  tape,  converted  to  an  analog 

XgjlaX,      cinQ     J- SQ      uO      ucimxiicixS      ±  H /.  u      <^  il  ^ '^  <~x^  1 1      j. -i  ^ .      Oii      T_ji.O      ciriclJLLJy 

board  (corresponding  to  channels  1  through  4).   The  signal 
could  then  be  monitored  on  various  recording  devices. 
Because  of  the  short  duration  of  the  signal  it  was  best 
displayed  on  the  available  pen  recorder  whose  input  terminals 
are  found  on  the  analog  board  (Figure  3) .   The  time  base 
of  the  signal  returned  from  the  tape  could  be  lengthened  or 
shortened  by  varying  the  value  of  the  input  parameter  NDFL . 
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